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Abstract 

A large class of technically non-chaotic systems, involving scatterings of light parti- 
cles by flat surfaces with sharp boundaries, is nonetheless characterized by complex 
random looking motion in phase space. For these systems one may define a gener- 
alized, Tsallis type dynamical entropy that increases linearly with time. It charac- 
terizes a maximal gain of information about the system that increases as a power 
of time. However, this entropy cannot be chosen independently from the choice of 
coarse graining lengths and it assigns positive dynamical entropies also to fully in- 
tegrable systems. By considering these dependencies in detail one usually will be 
able to distinguish weakly chaotic from fully integrable systems. 
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1 Introduction 



Dynamical entropies measure rates at whicli information may be gained by 
observing a system of interest subjected to some law of temporal evolution. For 
deterministic dynamics one considers the time evolution of an initial ensemble 
of systems distributed uniformly inside a small compact box of diameter e 
(typically a ball or a cube) within the phase space characterizing the possible 
states of the system. One then asks how many boxes of size e are at least 
needed to cover the image of the initial box after a time t has elapsed. If, 
in the limit e ^ 0, the natural logarithm of this number increases as htopt 
asymptotically for t ^ oo, the quantity htop is called the topological entropy. 
In physical applications usually the Kolmogorov- Sinai entropy, designated as 
Hks is more relevant. This is defined in similar way, but the covering set is 
not required to contain the full image of the initial box; it only has to cover 
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this up to a subset of measure fi{t) satisfying limt^ooA*(0 — 0. To interpret 
this, assume e is the minimal rcsohition with which one may distinguish phase 
space points by physical observation. Then by repeating the observations at 
regular time intervals the number of distinguishable initial boxes increases as 
exphxst for a typical large but not exhaustive number of observed systems 
and it increases as exp htopt for an exhaustive search in which also rare events 
are weighted properly for the full time range of observation. 

For chaotic dynamical systems that are closed (no points escape from phase 
space) Pesin's theorem equates the KS-entropy to the sum of the positive Lya- 
punov exponents. Intuitively this result is not hard to understand. Each posi- 
tive Lyapunov exponent Aj corresponds to an independent direction in phase 
space with a smooth expansion by a factor of roughly exp X^t over a time t. 
So the overall expansion factor, and thereby the number of boxes needed to 
cover the image of the initial box, increases as exp J2i \t- A priori one could 
expect though, that, at least in part, this may be offset by shrinkage in the 
contracting directions. However, this contraction typically takes place towards 
a fractal set with a striated structure on increasingly finer length scales. Once 
the typical length scale of separation between striation sheets falls below the 
box size e the number of boxes required to cover the set does not decrease 
any more. For comparison: consider the number of intervals of length e needed 
to cover successive approximations of the middle-third Cantor set. Once the 
length of the subintervals that are erased at the next iteration falls below e 
the number of covering intervals essentially does not change any more. The 
difference between topological and KS-entropy may be understood again by 
observing that on subsets of phase space of measure zero the Lyapunov expo- 
nents may differ from the "bulk" Lyapunov exponents, which for an ergodic 
system are constant almost everywhere. 

There exist dynamical systems that are not chaotic, in the ususal sense of 
having at least one positive Lyapunov exponent, but still show a fairly rapid 
decay to equilibrium, related to motion in phase space that appears very ran- 
dom. Simple examples of these are wind tree m,odels, in which light particles 
without any mutual interaction move elastically among fixed polygonal scat- 
terers. Lepri et al.[l] numerically studied such a model consisting of diamond 
shaped scatterers with corner angles that are irrational fractions of vr on a 
triangular lattice with bounded free paths between collisions. They found this 
model equilibrates rapidly on the projection to the unit cell in space times 
the unit sphere for the velocities, and on large time and length scales it is de- 
scribed by normal diffusion. Dettmann et al.[3] looked at Ehrenfests original 
wind tree model[2], in which the scatterers are randomly distributed parallel 
squares and the velocities of the light particles always point along the diag- 
onals of the squares. In this case the velocities equilibrate among the four 
allowed directions and, in finite volume the density of a cloud of light parti- 
cles becomes uniform. Again, on large time and length scales the dynamics 



2 



becomes diffusive, provided the scattcrcrs arc not allowed to overlap each 
other [4, 5, 6]. Another variation, studied by Li et al.[7], consists of a channel 
lined by triangular scatterers, through which light particles move, giving rise 
again to diffusive motion on large scales. Other non-chaotic systems leading 
to macroscopically stochastic looking motion are mixtures of point particles 
of different masses in one dimension. These, for particles, may be viewed as 
a type of wind tree model on a 2A^-dimensional phase space bounded by fiat 
hyperplanes. The intersections between these hyperplanes take the role of the 
corners; "hitting the other side of a corner" translates to an interchange in the 
sequence of two subsequent collisions. Livi et al.[8] and Grassberger et al.[9] 
found that such systems show anomalous hydrodynamics, which is typical for 
one-dimensional systems. As a last example I would like to mention systems 
of many moving parallel oriented polygons (in 2d) or polyhedra (in 3d) , which 
in coUisions exchange the velocity components normal to the plane of colli- 
sion. These systems too may be viewed as wind tree like models in a high 
dimensional phase space. No doubt they will show normal hydrodynamics (up 
to the usual nonlinearities in 2d, see e.g. [10]) although they are obviously 
non-chaotic. But to my knowledge no detailed investigations have been done 
for such systems so far. 

For the systems described above, according to Pesin's theorem the KS-entropy 
equals zero. Yet, it is clear that information may be gained by repeated ob- 
servation, only the number of distinguishable initial boxes will not increase 
exponentially with time, but rather as a power law. Zaslavsky and Edelman 
refer to this situation as pseudochaos[ll]. I will use the term "weak chaos" in 
the sequel, without giving a precise definition. For quantifying the power law 
increase of the information on the initial distribution Tsallis' formalism seems 
suitable, as it is precisely meant to define an entropy for systems in which 
the countable number of states (here distinguishable initial boxes) does not 
increase exponentially with system size (here time) but rather as a power law. 
But there are some subtleties here to deal with. As we will see, the dynami- 
cal entropy as defined usually, will become dependent on the coarse graining 
lengths in the system. It has to be multiplied by powers of these to allow for 
taking the limit e ^ 0. In addition, even very regular systems, such as an ideal 
gas in a periodic box, or a weakly anharmonic oscillator, will turn out to give 
rise to a positive dynamic Tsallis entropy, thereby excluding this property by 
itself as an indicator of weak chaos. Keeping track of the power law by which 
the information increases may help to recognize weak chaos. Even this requires 
a lot of care, because the powers of time resulting from integrable degrees of 
freedom may vary, depending on the type of motion they represent. 

In the next sections we will consider dynamical Tsallis entropies and look at 
wind tree models in detail. In the discussion the feasibility of using a dynamical 
TsaUis entropy to characterize weak chaos will be discussed further. 
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2 Dynamical Tsallis entropies 



Supose the phase space of a class of dynamical systems is divided into subsets 
Cj such that the probabihty of finding a system in Cj equals Pi. Then the 
TsaUis entropies are defined through. 

S, = ^ (1) 



Here g is a free parameter the choice of which should be dictated by the 
properties of the system considered. In the limit q ^ 1 the Tsallis entropy 
reduces to the usual entropy. 

In the case of dynamical systems evolving in time from the interior of an initial 
set of diameter e one may choose for the Cj a covering set of non-intersecting 
subsets of diameter e and define Pi{t) as the fraction of initial phase points 
that end up inside Cj after the time evolution over time t. The number of Cj 
with non- vanishing Pi{t) will increase as a power of t. Now let us suppose that 
on this subset pi{t) for large t behaves as 



with Wi distributed according to some probability distribution with a mean 
value of unity. Then, in order for Sq{t) to increase proportionally to t one has 
to choose q — 1 — 1/ a. In this case the entropy behaves as 

5i_i/„(i) = aCV"<^f-V")>, (3) 



with the brackets indicating an average over the distribution of Wi. In case 
this is an exponential distribution, p{w) = exp—w, one finds < >= 
r(2 - 1/a). 

A further remark is required here. As time increases one should choose increas- 
ingly smaller values of e in order to avoid Cj's containing parts of the image 
of several initial boxes. For fixed e this will always occur after sufficiently long 
time. In principle one should take a limit e — > followed by the limit t ^ oo. 
We will see however that the generalized dynamical entropies one would like to 
define are dependent on e, so one cannot take these limits in a straightforward 
way. In fact, as we will see, it can only be done after redefining the dynamical 
entropies in an appropriate way. 
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Fig. 1. A configuration of the wind tree model with randomly oriented square scat- 
terers. The dashed line represents a posible light-particle trajectory (note the peri- 
odic boundary conditions) 

3 Wind tree models 



Here I will consider a wind tree model in two dimensions with square scatterers 
of side length a that are randomly oriented and distributed randomly over a 
rectangle with periodic boundary conditions. The scatterers are not allowed 
to overlap each other. I will assume that the density p of scatterers is small, 
that is pa^ ^ 1, so subsequent collisions of a light particle are to a very 
good approximation uncorrelated. For the subsets Cj of the phase space let 
us choose small rectangular blocks of side lengths in position space and 
in velocity space. Over a time t the dynamics stretches out these blocks 
into parallellopipida. For large t, in the direction parallel to the velocity these 
obtain an extension of length e^t. In the direction orthogonal to the velocity 
the front width covered by the parallellopipidum also grows as Cyt. In addition 
the block is split into two subblocks, as illustrated in Fig. 2, as soon as its front 
hits a corner of one of the scatterers. Next, when one of the two subblocks 
hits another corner it is split up likewise and so on. The total front width of 
all the propagating subblocks increases with time as e^i, independent of the 
number of splittings that have taken place. Therefore the average number of 
splittings during a long time t is given by 



N{t) = -pvej\ 



(4) 



Here v is the average speed of a light particle. The average area swept out 
by the total front of all the subblocks is and this is multiplied by 

3p to obtain the average number of visible corners in this area. The average 
number of resulting subblocks is almost equal to N{t). To obtain the number 
of subsets C{i) needed to cover all of the subblocks at t one has to multiply 
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Fig. 2. When a block of nonzero density in phase space hits a corner of a scatterer 
it is spHt into two subblocks. If it just hits a side it is reflected in its entirety. 

N{t) by e^t/er, the number of subsets needed to cover the extension in the 
parallel direction. So in Eq. (2) a = 3 and C = 3pvel/2er- The positions of the 
corners with which the collisions take place follow a Poisson distribution along 
the front, hence the distribution of Wi is an exponential one indeed. Using the 
results of the preceding section we find 



4 Discussion 



The first thing that draws the attention in Eq. (5) is the dependence of the 
entropy on the coarse graining lengths e^, and e^-. This is a striking difference 
from the result found by Borges et al.[12] for the logistic map at the onset of 
chaos. Here too the number of boxes increases algebraically with time, but this 
is due to a power law separation between infinitesimally close trajectories in 
phase space. In the wind tree models infinitesimally close trajectories typically 
just separate linearly with time, due to the infinitesimal velocity difference, 
but the interesting power law behavior in the time dependence of the number 
of covering sets needed is caused by the corner induced splitting of trajectories 
at arbitrarily small, but always finite distances from each other. The same is 
true for all the other examples mentioned in the introduction. As mentioned 
at the end Sec. 2 a sharp definition of dynamical entropies involves taking a 
limit e ^ followed by the limit t ^ oo. Obviously this makes no sense for the 
entropy obtained in Eq. (5). The limit e — > is not well-defined. Apparently 
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the only way to repair this is by defining still another dynamical entropy as 



An alternative way of obtaining this is by requiring 6^ = 6^. For wind tree like 
models in d dimensions the prefactor e^/e^ should be raised to the power d—1. 
This solution is not very elegant, e.g. there seems to be no obvious choice for 
making the e's dimensionless. But at least it seems to give a proper limiting 
behavior for e — > for all wind tree like models. 

A second point of concern is the entropy of an ideal gas in a periodic box. One 
easily sees that, with the same type of covering sets as for the wind tree model, 
the number of sets needed to cover the time evolved initial set increases in d 
dimensions as [eyt/erY. So for g = 1 — l/d one now has a dynamical entropy 
that increases linearly with time. The same occurs for other integrable systems 
in which the periods of the angle variables depend on the values of the action 
variables. This is usually the case, with harmonic oscillators being the notable 
exception. For (i-dimensional wind tree models the power of t by which the 
number of needed covering sets increases with time is 2(i — 1. So it looks like 
one could use the condition q > 1 — 1/dasa criterion for weak chaos in a 
(i-dimensional hamiltonian system. Here too care is needed. For example in a 
three-dimensional system that is a wind tree model in the x — y plane but in 
addition performs a harmonic oscillation in the 2;-direction, one has g = 1/3, 
just as for a three-dimensional ideal gas. A distinction between the two cases 
may still be found in the dependence of S2/3 on and e^. The powers of 
and er by which one has to multiply Sg in order to obtain a S with a proper 
limit for e ^ obviously are different for integrable systems of various types, 
for wind tree like models and for mixed systems. 

It is not hard constructing generalizations of the wind tree model to higher 
dimensions. One just has to replace the scatterers by polyhedra of appropriate 
dimension. And indeed for each direction orthogonal to the velocity one picks 
up a factor t"^ in the growth factor of N{t). Similar observations can be made 
for the moving polygons. 

To conclude: Dynamical Tsallis entropies may be useful for characterizing 
weakly chaotic systems in which the information on initial conditions increases 
as a power law rather than an exponential. For wind tree like systems however, 
the precise value of this entropy depends on the choice of the coarse-graining 
lengths, which basically sets the smallest scales on which one may distinguish 
different systems in a single observation. To obtain a dynamical entropy with 
a proper limit when these coarse graining lengths are sent to zero, one first has 
to multiply the "straightforward" dynamical entropy with appropriate powers 
of the coarse graining lengths. Further, the generalized entropies typically 
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assume non-zero values also for integrable systems. In many cases the g- values 
needed to find an entropy that increases linearly with time will allow one to 
distinguish wind tree like systems from integrable ones, but this too is not 
always the case. If not, the dependence on the coarse graining lengths may be 
used for a further distinction. 

HvB acknowledges support by the Mathematical physics program of FOM and 
NWO/GBE. 
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